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We analyze the effects of quantum correlations, like entanglement and discord, on the efficiency 
of phase estimation by studying four quantum circuits that can be readily implemented using NMR 
techniques. These circuits define: a standard strategy of repeated single qubit measurements, a 
classical strategy where only classical correlations are allowed, and two quantum strategies where 
nonclassical correlations are allowed. In addition to counting space (number of qubits) and time 
(number of gates) requirements, we introduce mixedness as a key constraint of the experiment. We 
compare the efficiency of the four strategies as a function of the mixedness parameter. We find that 
the quantum strategy gives V^V enhancement over the standard strategy for the same amount of 
mixedness. This result applies even for highly-mixed states that have nonclassical correlations but 
no entanglement. 



I. INTRODUCTION 

There is a great deal of work on optimal phase estima- 
tion p] [2] addressing the practical problems of state gen- 
eration, particle loss and decoherence. However, this has 
mainly been done within specific experimental contexts 
and often with (initially) pure states of the probe only pQ . 
To understand the origin of the quantum enhancement 
over the standard quantum limit, many have analyzed 
the role of the number of bits required and the number of 
elementary gates needed, as well as the role of entangle- 
ment [SJ. However, in addition to counting the resources, 
constraints also need to be taken into account. For 
example, in nuclear-magnetic-resonance (NMR)-based 
quantum information processing, the quantum opera- 
tions take place at a fixed (room) temperature. This, 
of course, means that not all physical states can be ac- 
cessed, only those of a certain (fixed) degree of mixedness. 
When optimizing phase estimation, this mixedness has to 
be taken into account. In fact, the degree of mixedness 
now becomes at least as fundamental as the requirements 
of the number of qubits and gates. 

The other element that plays a crucial role is corre- 
lations, namely entanglement when dealing with pure 
states. However, quantifying correlations as a resource 
and mixedness as a constraint, leads to a complicated pic- 
ture. For mixed states entanglement is no longer the sole 
correlation present; other quantumness quantifiers like 
quantum discord [1HS] may be relevant. A well-studied 
example of this sort is the deterministic quantum com- 
putation with one qubit (DQC1) 7 . Here, a classically- 
hard task is performed efficiently quantum mechanically, 
but no (or only marginal) entanglement is present, while 
quantum discord can be present even when entanglement 
is vanishing; this led to the conjecture that discord maybe 
responsible for the quantum speed-up [5]. 



In this article, the role of correlations in quantum 
metrology is studied along the lines of [8]. We compare 
different strategies at a given (fixed) mixedness, within 
the constraint where pure states are not readily available 
and classical noise is always present (in contrast to the 
framework of [9J 110)). Our study is intended to gain in- 
sight into how mixed-state correlations, namely entangle- 
ment and discord, contribute to quantum enhancement. 
We show that mixed-state metrology leads to the same 
uncertainty in phase estimation as pure states but with 
an overhead that scales linearly with the classical noise. 
This turns out to be independent of entanglement, and 
therefore a quadratic quantum enhancement is available 
even for states that are highly mixed and fully separable. 



II. FRAMEWORK FOR CORRELATIONS 
STUDIES IN MIXED STATE METROLOGY 

We work with an iV-qubit system with each qubit ini- 
tially being in the mixed state 
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From this we construct correlated states (also called 
probe states) of various types with unitary gates. Re- 
call that global unitary operations preserve the mixed- 
ness of the total state but not the correlations con- 
tained within it. We study three strategies having dif- 
ferent types of multipartite correlations. The first two 
are quantum strategies, called Ql and Q2, which use 
GHZ-diagonal states. These states have quantum corre- 
lations such as entanglement and discord. The third is a 
classical strategy, labeled CI, which uses only classically- 
correlated states (defined as having zero discord) . We 
compare these three strategies to the standard strategy, 
called S, where a single qubit is used ./V times to estimate 
the phase. Below we lay out the details of preparing these 
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FIG. 1. The circuits for the four strategies considered in this 
paper are shown above, pi is the control qubit and pN-i = 
®il2 Pi are the rest of the N — 1 qubits. Initially all qubits 
are in the same state, given in Eq. [I] 



states. The circuits for preparing these states are explic- 
itly given in Figs. [TJa)-(d). 



A. States preparation 

1. Standard strategy 

The state for standard strategy is obtained by applying 
a Hadamard gate, H , to each qubit 



<8>AT 



qs = (HpH) 



2. Classical strategy 
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The classical state is created by applying a series of 
C-Not gates between the first and the ith qubit, Cu, 
followed by Hadamard gate on each qubit. 
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Above H N = pn-i = ®f =2 Pi and c = 

(S)il=2 where Cu is a C-Not operation with first qubit 
as the control and ith qubit as the target. 



3. Quantum strategy 1 

For the first quantum strategy the GHZ diagonal state 
is prepared by taking the initial uncorrelated TV qubit 



state and applying the Hadamard gate to the first qubit 
followed by a series of C-Not gates between the first and 
the ith qubit. 



qqi =CHt [pi ® pn-i) HiC 

_1 / p®*- 1 pipa^' 1 
-2 [vfap)*"- 1 (axpa.)®"- 1 



(4) 



Above Hi = H ( 

This state was employed in the experiment reported 
in but it turns out not to be the optimal state. A 
more optimal state is described below. 



4- Quantum strategy 2 

For the second quantum strategy the GHZ diagonal 
state is prepared by taking the initial uncorrelated N 
qubit state and applying C-Not gates between the first 
and the ith qubit followed by the Hadamard gate to the 
first qubit followed by another series of C-Not gates be- 
tween the first and the ith qubit. 

qq 2 =CHxC (pi (g) pn-i) C HiC 

_W p^' 1 (pa x r N ^ 
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The Q2 state is constructed in much of the same way 
as the Ql state, but initialized with C-Not gates to shift 
the initial population. This strategy was used in the 
experiment reported in [12] . 



B. Relations to NMR 

Our model corresponds particularly well to a set of re- 
cent NMR experiments on magnetic-field sensing [TT1I12) . 
The initial state in NMR experiments is 'pseudo-pure' — 
a density matrix which is very close to being completely 
mixed although its eigenvalues are not quite identical. 
In NMR, the qubits are the spins of nuclei and unitary 
operations on these spins are performed by applying elec- 
tromagnetic pulses of a selected frequency and duration. 
Qubits can be selectively addressed by choosing spins 
with specific resonance frequency; these can be local or 
global (entangling) unitary operations. As more species 
of spin are added the pulses needed to exclusively address 
and couple the different species becomes more difficult. 
However in practice, these operations can be performed 
with extremely high fidelity (see [13] for detailed analy- 
sis). 

In the experiments reported in [TTJ [T^] only two species 
were used. This is the so called star topology, where the 
first qubit (pi in Fig. ffl) is used as the control qubit and 
the rest (pn-i m Fig.^J are subjected to a single trans- 
formation at once. For us this translates into using the 
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same one-qubit gate on each of the qubits in pn-i and 
a two-qubit gate between the control qubit and each of 
the qubits in Pn-i- The state preparation in [TT] slightly 
differs from the state preparation in |12) . The difference 
is precisely the difference in the two quantum strategies 
considered here: the state in [TT] corresponds to the cir- 
cuit in Fig. [ljc) and the state in [T2] corresponds to the 
circuit in Fig. [ijd). 



III. QUANTUM FISHER INFORMATION FOR 
DIFFERENT STRATEGIES 

Now we are in the position to compute the phase un- 
certainty for each the strategy above. For mixed states, 
the phase uncertainty is determined by computing the 
quantum Fisher information [T4HT6] given by 
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where {rjj} and {l^j)} are the eigenvalues and the cor- 
responding eigenvectors of state g, and G is the Hamil- 
tonian of the process that the state is subjected to. The 
Hamiltonian for the ith party is Gi = For the 

iV-party case, each party picks up a phase locally, which 
means that the global Hamiltonian is G = Gi <E) 1^, 
where the identity matrix acts on the remainder of the 
Hilbert space. The phase uncertainty is related to quan- 
tum Fisher information as 



This is the same as doing the phase estimation experi- 
ment with a single qubit and repeating the experiment N 
times. The initial state of the qubit is taken to be g$ — 
(HpH)® N . The eigenvectors are = |n, . . . , rjv) , 

where |r») € {|+) , |— )} is the eigenstate of the ith sub- 
system. We denote an arbitrary degenerate eigenvector, 
having ( )-fold degeneracy, as \ip m ): Label m counts the 
number of subsystems in state |— ). The corresponding 
eigenvalue is x] m — 

Now, we label the eigenvectors of qubits 2 to N by 
|x) and consider the eigenvectors \ip m ) = \+,Xm) and 
\tp m +i) = | — ,Xm) and the action of the Hamiltonian on 
them (ip m \G\ip m+ i) = Y!,i(+>Xm\Gi\-,Xm)- The only 
term that survives is (+|Gi|— ) (Xm\Xm) — §■ Since the 
states are in the product form, the same result is true for 
all subsystems and the quantum Fisher information for 
N qubits is N times the quantum Fisher information of 
a single qubit 
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This is the expected result and agrees with the pure state 
results asp-> 1. 
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Quantum Fisher information is a function of the 
Hamiltonian that generates the interaction between the 
probe and the object being measured. It also depends on 
the state of the probe. In our problem, the Hamiltonian 
is the same for all strategies, only the correlations within 
the states change. The final measurements at the end 
are assumed to be optimal generalized measurements as 
is assumed in the derivation of the quantum Fisher infor- 
mation. For the strategies considered here the measure- 
ments turn out to be rather straight forward, see |13j for 
details. The equality in the last equation can be achieved 
by statistical estimators provided the system is sampled 
several times. A detailed analysis would identify a statis- 
tical estimator to extract the maximum information and 
saturate the Cramer- Rao bound [17 . We compute the 
quantum Fisher information and the phase uncertainties 
for the three strategies discussed here as follows. At the 
end of the section we compare these values. 



A. Standard strategy 

We begin by computing the quantum Fisher informa- 
tion for N qubits that share no correlations whatsoever. 



B. Classical strategy 

To create a classical state we start with p® N , which 
has eigenvectors \^j) — \r\, . . . , rjy) , where \n) — |0) , |1) 
is the eigenstate of the ith subsystems. Once again we 
denote an arbitrary degenerate eigenvector, having ( V 
fold degeneracy, as \tp m )' Label m counts the number of 
subsystems in state 1 1) . The corresponding eigenvalue is 

Vm — A Q A i ■ 

Next we apply the C-Not gate between the first and 
the ith qubit. The eigenstates under the C-Not oper- 
ation change as following: \ip m ) = \0,Xm) -> \0,Xm) 
and |VWfi) = \hXrn) \ipN-m) = |1, XJV-m-i)- Next 
a Hadamard gate is applied on each qubit, which sim- 
ply changes |0) — > |+) and |1) — > |— ). The action of 
the Hamiltonian on the eigenstates \ + ,Xm) and | — ,Xm) 
gives (+ , Xm\G\ — , Xm) = \ with the corresponding left 
and right eigenvalues r]i = A^ _m Af and r) r = A™Af ~"\ 



\Xm) occur with binomial distribution ( N m 1 ) 

Hamiltonian 



of the 
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is 



The action 
eigenstates \±,+i,x m ,i 

(±,+i,Xm,i\G\±,-i,Xm,i) = (+i\Gi\-i) = \ with 
the ith state on the left is |+) and the state on the right 
is I - )■ \Xm,i) is the state of parties excluding the first 
and the ith qubits occurring ( N ~ 2 ) times. The index i 
runs up to N — 1 yielding the same inner-product. 
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When the first qubit is in state |+) the correspond- 
ing left and right eigenvalues are r\i — A^~ m A" 1 and 
Vr = A _m_1 A5™ +1 . The difference of these eigenvalues 
squared divided by their sum is simply A^~ m-1 A™p 2 . 
When the first qubit is in state |— ) the correspond- 
ing left and right eigenvalues are r\i — A^A^ - ™ and 
rj r = AQ i+1 Af The difference of these eigenvalues 

squared divided by their sum is simply A™ A 1 ^~ m ~ 1 p 2 . 

Due to symmetry all other Hamiltonians will have the 
same result as above. The Fisher information is simply 
the sum of the three results above 
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^Np 2 + l-p 2 -e- Np \ (9) 
where the last approximation is obtained numerically. 

C. Quantum strategy 1 

The eigenstates of qq\ are of the form \ip m ±) = 
^ (|0,Xm) ± |1, o"® JV_1 Xm)) with the + eigenvalues be- 



ing A^ m A™ and the — eigenvalues being Aq 
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and m denotes the number of subsystems in state 
Once again the degeneracies follow the binomial distri- 
bution. The action of the Hamiltonian is 

(V> m -|G|VW) =^«0,Xm| + (l,af N - 1 Xm \) 

xG(\0 )Xm )-\l,of N - 1 x m )) 
=i((0, Xm | + <l,af JV - 1 Xm |) 

x (m|0,x m }-(iV-m)|l,^ w - 1 x m » 



=-(2m-N). 
The quantum Fisher information is 
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TABLE I. Quantum Fisher information. The quantum Fisher 
information gives the lower bound on the phase uncertainty, 
> 1/ \/F(q). Note that classical noise in strategies £ and 
are the same, p 2 , independent of TV. The phase uncer- 
tainty for each strategy is plotted as a function of p in Fig. [2] 



Once again, the result above satisfies the known result 
for pure states. Note that the leading term goes as p 2 N 
and N 2 term has a pre- factor of p A . 



D. Quantum strategy 2 

The eigenstates of qq2 are of the form \tjj m ±) = 
(|0,Xm> ± l^of^Xm)) with the + eigenvalues be- 
ing A Y_m A™ and the - eigenvalues being Ag l Af r_m and 
m = 0, . . . , TV denotes the number of subsystems in state 
|1). Note that the eigenstates here are the same as the 
previous case but the corresponding eigenvalues are dif- 
ferent. Once again the degeneracies follow the binomial 
distribution. The action of the Hamiltonian is same as 
the previous case. The quantum Fisher information is 
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Once again, the result above satisfies the known result for 
pure states. Numerical results indicate that F(qq 2 ) > 
p 2 N 2 . This means that the classical noise is the same as 
the standard case but we have a quadratic enhancement 
in the number of qubits. 



E. Comparison of quantum Fisher information of 
different strategies 

The results of phase uncertainties, presented in the 
Table |IJ are plotted in Fig. [2] for TV = 10. Remember 
that our goal is to compare different strategies at a fixed 
mixedness, i.e. a fixed value of p, while changing the 
number of qubits, i.e. the value of TV, does not change 
the overall behavior of these curves. Strategy Q2 is far 
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rigorous proof of the optimization of the quantum Fisher 
information for all N and p is a very difficult problem. 
Below we argue that for any N the states chosen for the 
strategies CI and Q2 are optimal for p close to 1 and we 
conjecture that they remain so for all values of p. Cer- 
tainly they provide strong lower bounds for the quantum 
Fisher information sufficient to support the conclusions 
of this article. We should reemphasize these probe states 
are experimentally realizable and realistic. 



A. Optimizing the standard strategy 



FIG. 2. (Color online.) Phase uncertainties for N — 10. 
The lower bound on the phase uncertainty as given by the 
quantum Fisher information as functions of p (see Table |T|. 
Q2 is by far the most optimal strategy for all values of p. 
While, CI strategy is better than Ql for small values of p. 



better than any of the other strategies, especially Ql. 
In fact, for highly- mixed states the CI strategy is better 
than Ql. The point at which the classical strategy over- 
takes Ql strategy is approximately whenp ss (This 
crossing point turns out to be independent of entangle- 
ment as the crossing occurs before (for N = 2) and after 
(for N > 2) entanglement vanishes.) Quantum Fisher 
information for strategy Q2 is affected by the classical 
noise in qubits only quadratically, i.e. F(gQ 2 ) ~ (Np) 2 ; 
it could have been exponential in the number of qubits, 
i.e. F(qq 2 ) ~ N 2 p N . Photon losses for optical setups 
have a devastating effect on quantum enhancement, while 
the NMR setup seems to be robust under lack of initial 
coherence. 



For the standard strategy, the quantum Fisher infor- 
mation can be computed for a single party and the N- 
party quantum Fisher information is simply N times the 
former. The single-party Hamiltonian for the process is 
G = |1) (1|. Therefore the eigenbasis for the density ma- 
trix should be {|+) , |—)} to maximize the transition from 
one eigenstate to another. This is why the Hadamard 
gate is applied on all qubits for the preparation. 



B. Optimizing the classical strategy 

A classical state has a separable (locally orthonormal) 
eigenbasis [5], therefore a classical state is simply ob- 
tained by rearranging the correspondence between the 
eigenvectors and the eigenvalues of the 7V-qubit density 
matrix of the standard strategy (Eq. [2]) . Therefore the 
unitary operations can only permute the computational 
basis along with local rotations. 

The eigenvectors of the classical state are given as 



m) = ri r 2 . 



■ r N ) 



The action of the Hamiltonian is 



IV. OPTIMALITY AND BOUNDS 

The quantum Fisher information in Eq. [6] is a func- 
tion of how the process Hamiltonian can connect two 
eigenstates of the density matrix and the difference in 
the corresponding eigenvalues. Maximizing the two will 
maximize the quantum Fisher information subject to the 
constraint of the correlation class. Since only unitary op- 
erations are allowed for preparation, the spectrum of the 
density operator remains hxed for all strategies. There- 
fore the first term of Eq. [6j i.e. (jij — r/k) 2 '/(f]j + Vk) is 
fixed. The only change can come from the changes in the 
eigenvectors. The optimal quantum Fisher information 
is then given by 



F maK (g) = 4max^ 

3> 



(13) 



where the unitary transformation U has to be constrained 
such that it does not change the correlation class of g. A 



G^j) =Y,Gi\ri r 2 ...r N ) 

i 

= Y,(Mn)\n...U...r N ). (14) 

i 

For (vpfclGl^j) to be nonvanishing, |^) must only differ 
from only at one site. Then (* fe |G|<F,-) = | (li\n) | 2 , 
and the maximum is attained when |r) G Since 
the process Hamiltonian is diagonal in z— basis, we would 
like to rotate the eigenvectors to the x— basis by apply- 
ing the Hadamard gate to each qubit, similarly to the 
standard strategy above. 

Now we provide a prescription for maximizing the 
quantum Fisher information, which is certainly optimal 
for large p, and we conjecture that it remains optimal 
for all values. The key idea is to insure that the largest 
and smallest eigenvalues are connected by the process 
Hamiltonian, i.e. 



max - 



(Vj -Vk) 2 
Vj + Vk 



V j,k 



(15) 
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with (y>j\G\^ k ) ^ 0. The largest eigenvalue is A^, be- 
longing to eigenvector !+•••+)■ The action of the Hamil- 
tonian is on 1+ • • • +) is 



G\ 



N 



71 



(16) 



This is a superposition of N terms with the ith party 
in state |1). Therefore the only states that have a finite 
value for (ip m \G\ + ...+) ^ are: (1) the state with the 
lowest eigenvalue (Af), which becomes 



C-Not: |- -...-)-> | 
(2) states with one excitation, i.e. 



(17) 



. . . +) (there 

are N — 1 such states with eigenvalues A^ _1 Ai). These 
latter contributions will be small because the eigenvalues 
will be different by only one excitation, but occur mul- 
tiple times. Therefore, A^ — Af is the largest possible 
leading term for the quantum Fisher information. The 
same argument can be repeated for the eigenvectors with 
the second largest and the second smallest eigenvalues, 
and so on, until all eigenvectors are paired in this manner. 



Let us now show that the last term is the largest pos- 
sible leading term. Suppose that the eigenvector with 
the largest eigenvalue is connected to some other state 
not having the smallest eigenvalue and this is the leading 
term. Explicitly, we have 



MJV \N— m\m\2 
[A Q — A 1 A Q ) 

A jv + A jv-r» A „ t I (VtaaxIGIV* 
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Since the last term in the last equation is independent of 
p, when we take p = lwe have 



(Vtaax|G|V> m ) | 2 < ^V 2 
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N 2 is the largest possible value for quantum Fisher in- 
formation. We also have 



(A^-Af—A^) 2 ^ (Xq - Af ) 
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with equality satisfied if and only if m = when p ^ 
0, 1. This is true because the numerator becomes smaller 
and the denominator becomes larger as the value of m 
increase. Therefore 



C. Optimizing the quantum strategy 

We know for case p = 1 the optimal pure quantum 
state for metrology is the GHZ state. For the case p < 1 
we conjecture that a GHZ basis is the optimal basis for 
quantum Fisher information. Our first attempt along 
this lines is to use the same circuit as would be used 
for the pure state case, transforming the eigenstate with 
the largest eigenvalue, A^, into the GHZ state. This in 
fact strategy Ql. The problem with this strategy is that 
the process Hamiltonian connects this state to a second 
state that has an eigenvalue that is different by only one 
excitation, i.e. A^ _1 Ai. More precisely, the first term of 
the quantum Fisher information for strategy Ql is 



(A^ - Af r - m Ag 1 ) 2 

+ Af" m A; 



W-A 



JV-l 



A '>V 
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However, following line of reasoning of the Sec. IV B 
we would like the two eigenvalues for the leading term 
to be maximum and minimum. Therefore, it would be 
desirable to permute the eigenvalues of all eigenstates 
whose leading term is |1). This is precisely what the 
initial C-Not gates do in Q2: 

CHiC |11 ... 1) =CHi 1 10 ... 0) = -pC |-0 .. . 0) 

v 2 



= - 7 =(|00...0)-|ll...l)) (19) 

The leading term of quantum Fisher information for 
strategy Q2 is then 



(A w -Af) 2 2 
A~ + Af ■ 



(20) 
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(A "-Af) 2 



\q + A 



N- 



■N 1 



(24) 



with equality for m = 0. 

Now that the largest and the smallest eigenvalues are 
taken care of, we can repeat this process, matching the 
mth smallest eigenvalue with mth largest eigenvalue. 
These arguments strongly suggest that the strategy Q2 
is the optimal quantum strategy for mixed states for all 
values of p and N. 



V. CORRELATIONS FOR DIFFERENT 
STRATEGIES 

In order to relate the results of phase estimation to cor- 
relations, we have computed the all of the correlations for 
the strategies C/, Ql, and Q2. The state in the strategy S 
has no correlations and the state in CI strategy does not 
have any entanglement or discord by definition. Strate- 
gies Ql and Q2 have entanglement for some values of p, 
while quantum discord is present for all values of p. We 
begin with computing the entanglement vanishing points 
for Ql and Q2. 



A. Entanglement in qqi 

It has been shown that a necessary and sufficient con- 
dition for a GHZ-diagonal state to be separable is that 
every possible partial transposition is positive [18] . Us- 
ing this result we can find a relation for a given iV that 
gives the value p for the boundary between separable 
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and entangled. The form of the states we are looking 
at in the computational basis have already been given in 
Eq. |4j Since this matrix is a collection of 2 x 2 block 
matrices, the partial transposition that gives the most 
negative eigenvalues is simply the one that results in a 
2x2 matrix with the smallest diagonal elements and the 
largest off-diagonal elements. Assuming that p > the 
2x2 matrix with the smallest diagonal elements is the 
one spanning the space in the central part of the matrix 
with diagonal elements A^" 1 . The 2x2 matrix with 
the largest off-diagonal elements sits in the corners with 
values pXg^ 1 . There exists a partial transposition that 
swaps the smallest off-diagonal elements with the largest 
off-diagonal elements resulting in the matrix 



N=5 



Aj P An 



p A 



N-l 



Af" 



(25) 



The smallest eigenvalue is then (A^ -1 — pA^ _1 ) and it 
is zero (this is the point the state becomes separable) at 
A^ _1 = pXq^ 1 . One can solve this equation numerically 
for a given TV. 




B. Entanglement in QQ2 

Using the same technique as above and assuming that 
p > in Eq.[5]we can show that the 2x2 matrix with the 
smallest diagonal elements has the elements IAq'A-l" 1 ", 
for N ± = N/2 for even N and N ± = (N±l)/2 for odd N. 
The 2x2 matrix with the largest off-diagonal elements 
sits in the corners with values A^ — Xf . There exists a 
partial transposition that swaps the smallest off-diagonal 
elements with the largest off-diagonal elements resulting 
in the matrix 



2A^-Af+ 



. N- 



Af 



\N \N o\ iv -\ Ar + 
A - A 1 ZA Q A-l 



(26) 



The smallest eigenvalue is then ^2A^~ A^ + — A^ + A^ 
which is zero (this is the point the state becomes sepa- 
rable) at A^A^ + = (A^ — A^) /2. Again, one can solve 
this equation numerically for a given N. 



C. Discord in qqi 

Quantum discord, denoted by D, is defined as the dis- 
tance (using relative entropy) between a quantum state 
and it's closest classical state: D(g) = miii{|k)} >S , ((?||x e ). 
The closest classical state, x e , is found by dephasing 
g in a locally-orthonormal product basis {|k)}: Xg = 
X; k |k)(k|e|k)(k|, (see [5] for details). We note that 
discord serves as the upper bound on entanglement as a 
function of p, E(p) < D(p) [6]. 

Computing quantum discord is an extremely hard 
problem; there exists no closed form solution even for 



FIG. 3. (Color online.) Numerical simulation of discord for 
Q1(Q2). Simulation of quantum discord (D) for TV = 5 with 
10,000 random measurements on Pqi(q2) to obtain Xqi(Q2)- 
The top (blue) line is the maximum possible values for D, 
i.e. 5 — S(qqi(Q2))- The bottom (red) line is the conjectured 
D = S(xqi(Q2)) - S{qq 1{Q2 ))- All random measurements 
points fall between the two lines suggesting that the conjec- 
tured formulae to be correct given Table [IT] 



arbitrary two-qubit states: The main difficulty lies in de- 
termining the minimizing basis {|k)}. In this problem 
we are dealing with a multi-qubit state. Using relative 
entropy of discord avoids making arbitrary bipartitions 
as would be required for computing bipartite discord. 

However, following the recipe of [B] , the closest classical 
state to qqi is conjectured to be given by dephasing qqi 
in the standard basis: 



XQi 



1 Ip^K- 1 

2 I 







(27) 



To calculate D we just need to take the difference in the 
entropies of qqi and XQi- 



D Qi = S(XQi) ~ S(g Q i) = 1 - S(p), 



(28) 



where S is the von Neumann entropy: S(p) — 
— tr[plog(jo)] and p is the state given in Eq. [T] 

Since the last equation is a conjecture, we have nu- 
merically simulated the closest classical state for up to 
five qubits (see Fig. [3]). The result above holds up (i.e. 
discord is independent of N), but we do not yet have 
an analytic proof. One can consider this result to be at 



Entanglement vanishing points 



E Q1 = O P<(\i/\o) N ~' 
E Q2 =0 o A^" Af + > 
Quantum discord 



Af 



Dqi 
D 



1 



S(P) 



Q2 - 2 E m=0 (m-l) ft 



Classical correlations 



A^+A^A. 



NS(p) 



Cat = (N - 1) [ft (Ag + A?) + ft (2A Ai) - S(p)] 
<7 Q1 = (JV-l)[l-S(p)] 
C Q2 = N[l - S(p)] - D Q2 

TABLE II. Above N± = (iV±l)/2 if AT is odd and N± = N/2 
when A^ is even. S(p) = — tr[plog(p)] is the von Neumann 
entropy and h(x) — —xlog(x). The expressions for discord 
(and therefore classical correlations) above are conjectured 
and supported by numerical studies (see Fig. |3j. The details 
of the calculations are below. We have plotted these correla- 
tions as functions of p in Figs. 4(a) and 4(b)| 



least an upper bound on discord. We should note that 
the lower bound discord is strictly greater than 0, as it 
is easy to verify it is a quantum correlated state [TS]. Fi- 
nally, we have plotted discord given in the last equation 
as a function of p in Fig. 4(a) 
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D. Discord in qq2 

Since both qqi and qqi are GHZ-diagonal states, the 
form of their closest classical states are also the same. 
Which means we can simply dephase QQ2 in the compu- 
tation basis to get: 



XQ2 



A / o^- 1 







V (a x pa x r N -\ 
Ai / Kpa,)^- 1 
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(8»JV— 1 



(29) 



To calculate D we just need to take the difference in the 
entropies of Qghz and X- 

Dq2 =S{xq 2 ) - S(q Q2 ) 

( \N— m\m _i_ \m\N—m\ 

=2 ^2h(^ Ai + A ° Ai 



X 1 m j - NS ^' 



(30) 



where h(x) = —a; log (a:). Once again this formula is con- 
jectured, but numerical evidence shown in Fig.|3]supports 
this result. Finally, we have plotted discord given in the 
last equation as a function of p in Fig. 4(a)| 



FIG. 4. (Color online.) Conjectured discord and classical cor- 
relations as functions of p. (a) Discord is always present for 
the two quantum strategies. Discord in Ql is independent of 
N. Entanglement is always (equal at p — 1 or) smaller than 
discord and vanishes (for = 10) around p = 0.118 for Ql 
and around p = 0.088 for Q2. (b) Classical correlations for 
N = 10 are plotted as function of the mixedness. 



E. Review of correlations in probe states. 

Finally, we compute classical correlations following 
again the procedure in [6] and shown in Table [II] We 
have plotted the classical correlations given in the last 
equations as a function of p in Fig. [4(b)] In Table [IT] 
we list the formulae for all correlations computed in this 
section. 



VI. ANALYSIS 

Now we are in the position to relate correlations with 
the enhancement of the quantum Fisher information. We 
start by noticing that quantum Fisher information is af- 
fected by classical noise in qubits only quadratically, i.e. 
F{QQ2) ~ {Np) 2 , while entanglement for Q2 vanishes 
when p ss 1 /N. Classical correlations for the three strate- 
gies scale linearly with the number of qubits N (see Ta- 
ble. ITTl). In fact, Ql has more classical correlations than 
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CI and Q2 for all values of p. This supports the expected 
result that classical correlations, although, present in 
bulk do not contribute to quantum enhancement. The 
total correlations, defined as the sum of quantum discord 
and classical correlations: T = D + C, are the same for 
both Ql and Q2. This further allows us to distinguishes 
the role of quantum correlations in the two cases. 

Which brings us to our main observations. The en- 
hancement of phase uncertainty (hence quantum Fisher 
information) due to the optimal quantum strategy over 
the standard strategy is 



Quantum Advantage = \l ~ VN V p. (31) 
V 

Since the classical noise is roughly p 2 for both strategies 
S and Q2 (see Table 0, the quantum advantage is VJV. 
This is true for highly-mixed states that have no entan- 
glement, i.e. p close to zero. Surprisingly, not a great 
deal of quantum coherence is needed to attain quantum 
advantage in quantum metrology. 

For the experiments reported in [Til El P ~ 10~ 5 and 
N ss 10. Therefore both states qq\ or qq 2 are unen- 
tangled. Both experiments reported quantum enhance- 
ment, which is in accordance with our findings. Quan- 
tum discord, on the other hand, does not vanish until 
p — > 0. And for Q2, quantum discord depends on the 



number of qubits, unlike for Ql (see Table. |Tl| ), assum- 
ing our conjectured expressions for discord are fully valid. 
Quantum discord for Q2 grows for small values of p as 
N increases. This provides evidence that quantum dis- 
cord may have some responsibility for the enhancement 
in quantum metrology. We should note that for entan- 
glement to appear when p ss 10~ 5 , the number of qubits 
has to be roughly N ss 10 5 . 

In conclusion we have analyzed the role of quantum 
and classical correlations in mixed-state phase estima- 
tion. We found evidence that classical correlations do 
not play a large role in quantum enhancement, as ex- 
pected. However, we also showed that quadratic quan- 
tum enhancement does not vanish as entanglement van- 
ishes. For such states quantum discord is present and is 
a growing function of the number of qubits. This adds to 
the evidence that quantum discord may be responsible 
for some quantum enhancements. 
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